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REAL-GAS EFFECTS ON THE DRAG AND TRAJECTORIES 
OF A NONUFTING 140° CONICAL AEROSHELL 
DURING MARS ENTRY 

By Dennis O. Allison and Percy J. Bobbitt 
Langley Research Center 

SUMMARY 

The present paper is a study of the influence of real-gas effects on the drag and 
trajectories of a 140° conical aeroshell durii^ entry into five model Mars atmospheres. 
Drag coefficients are obtained by using an approximate technique which enables their 
determination at every point along a trajectory with reasonable computer time. The 
standard point mass trajectory equations are used and a spectrum of entry conditions 
are considered that consist of velocities of 4.27, 4.57, and 4.88 km/sec (14 000, 15 000, 
and 16 000 ft/sec) and entry flight -path angles of 13.5°, 15°, and 17.5°. All calculations 
were made for an entry spacecraft with a ratio of mass to maximum cross-sectional area 
of 76.5 feg/m2 (0.487 slug/ft^). It is shown that a single curve may be chosen to repre- 
sent the drag coefficient variations within 2 percent accuracy for all five model atmo- 
spheres and all the entry conditions considered. 

INTRODUCTION 

Most of the mission analysis research on the atmospheric-entry phase of a Mars 
lander mission has been done with a constant ballistic parameter. (For example, see 
refs. 1 and 2.) Since no mass losses or area chaises are assumed in these analyses, a 
constant ballistic parameter implies a constant drag coefficient. Ordinarily, this constant 
drag coefficient is taken to be the hypersonic value measured in wind tunnels with air as 
a test medium. In reality, however, the drag coefficient of a body changes during entry 
and the manner in which it changes depends on the vehicle speed and the atmospheric com- 
position, pressure, and temperature. These drag changes result from the excitation of 
the higher degrees of freedom of the gas, such as vibration and dissociation (commonly 
called real -gas effects), in the shock layer ahead of the body where the gas is much 
denser and hotter than in the ambient stream. For a predominantly CO 2 gas mixture 
such as the Mars atmosphere, real-gas effects on the drag are larger than those encoun- 
tered in air for the same velocity and ambient conditions; they can have a significant 
effect on an entry vehicle's trajectory. 





The purpose of the present paper is to assess the magnitude of the real-gas effects 
on the drag coefficient for a 140° conical aeroshell during entry into five model Mars 
atmospheres.* The model Mars atmospheres are designated and defined as follows: 


Max Hp atmosphere 
Max pg atmosphere 
Mean atmosphere 
Min Hp atmosphere 
Min Pg atmosphere 


atmosphere with maximum average density scale height 
atmosphere with maximum density near surface 
mean model atmosphere 

atmosphere with minimum average density scale height 
atmosphere with minimum density near surface 


Drag coefficients are obtained by using an approximate technique which enables their 
determination at every point along a trajectory with reasonable computer time. The 
standard point mass trajectory equations (see appendix A) are used and a spectrum of 
orbital entry conditions are considered at an altitude of 244 km (800 000 ft) including 
velocities of 4.27 to 4.88 km/sec (14 000 to 16 000 ft/sec) and entry flight -path angles of 
13.5° to 17.5°. The entry spacecraft is assumed to have a ratio of mass to maximum 
cross-sectional area of 76.5 kg/m2 (0.487 slug/ft2) in all calculations. 


SYMBOLS 


Trajectory parameters were originally chosen in the U.S. Customary Units and the 
model atmospheres selected were in the SI Units. However, these quantities are expressed 
in both systems of units for the convenience of the reader. 

A maximum cross-sectional area of the aeroshell, m^ (ft2) 

_ , • 4 . Drag force 

Ct) drag coefficient, -= — - — - — 

ip V 2 a 

g acceleration of gravity, m/sec2 (ft/sec2) 

gg acceleration of gravity near the surface of Mars, m/sec2 (ft/sec2) 


Center 


Mars model atmospheres developed by George P. Wood of Langley Research 
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hoo free-stream static enthalpy per unit mass, m^/sec^ (ft^/sec^) 

hg static enthalpy per unit mass behind the shock, rtfi/sec^ (ft^/sec^) 

H altitude, km (ft) 

M mass of entry spacecraft, kg (slugs) 

Moo free-stream Mach number 

Mu molecular mass of undissociated gas mixture 

Nco 2 number of moles of carbon dioxide per unit mass in undissociated gas 

mixture 

Nj^g Reynolds number based on base diameter of model 

p_^ free-stream pressure, N/m^ (lb/ft2) 

p^ reference pressure, 0.1013 MN/m^ (2116 Ib/ft^) 

p_ pressure behind shock, N/m^ (Ib/ft^) 

free-stream dynamic pressure , hpoo^oo^ > N/m^ (Ib/ft^) 

a 

Tg mean radius of Mars, km (ft) 

Rq universal gas constant, 8315 kg-mVsec2-*^K-(kg-mole) 

(49 720 lb-ftVsec2-°R- (lb-mole)) 

Ru gas constant for undissociated gas mixture, Rq/Mu, m^/sec^-^^K 

(ft2/sec2-°R) 

t time, sec 

Too free-stream temperature, (®R) 

Tq reference temperature, 273.2*^ K (491.7® R) 
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Tg temperature behind shock, °K (°R) 

free-stream velocity component normal to shock wave, m/sec (ft/sec) 

Ug normal component of velocity behind shock, m/sec (ft/sec) 

Voo free-stream velocity, m/sec (ft/sec) 

Vg entry velocity, m/sec (ft/sec) 

y ratio of specific heats 

yg effective value of y, defined by equation (1) 

At time increment, sec 

e shock inclination angle, deg 

free-stream density, kg/m^ (slugs/ft 
pg reference density, MuP^^RoTo, kg/m^ (slugs/ft^) 

Pg density behind shock, kg/m^ (slugs/ft^) 

4> flight-path angle, deg 

entry flight-path angle, deg 

DRAG COEFFICIENT DETERMINATION 

A value of the drag coefficient Cp is required in the equations of motion (appen- 
dix A) in order to determine the position and velocity of an entry body after an increment 
of time At. Figure 1 presents experimental drag coefficient data (refs. 3 and 4) for 
Mach numbers from 2 to 10 for a 140° conical aeroshell obtained in two air facilities 
(y = 1.4). It shows that the drag coefficient is essentially constant and implies that there 
are no significant real-gas effects in air up to a Mach number of 10. However, in a 
Mars-like atmosphere composed largely of carbon dioxide, real-gas effects are e 2 q>ected 
to be important even in this speed range. In the remainder of this section, a procedure is 
outlined whereby approximate real-gas corrections to air-drag data may be applied for 
various atmospheric compositions. 
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Figtire 1.- Experimental variation of drag coefficient with 
Mach number for a l4o° conical aeroshell in air. 

There are three basic steps in the present approach for obtaining a real-gas drag 
coefficient for the model Mars atmosphere being "entered." The first step uses the 
instantaneous velocity, free-stream conditions, and equilibrium thermodynamic proper- 
ties behind the shock in solving the oblique-shock equations to obtain the real-gas shock 
density ratio. (See appendix B.) The second step uses this shock density ratio, a pre- 
scribed shock inclination angle (taken to be 70°), and free-stream Mach number to deter- 
mine an effective value of the ratio of specific heats from the perfect-gas relation 



(Note that eq. (1) is the solution obtained from the equation given in ref. 5 for the 
perfect-gas density jump across an oblique shock.) In the third step, the value of is 
used to determine the drag coefficient from a curve provided in figure 2. 

The variation of Cjj with y^ in figure 2 is based on a theoretical perfect -gas 
drag curve obtained by using the method of reference 6 for a 140° conical aeroshell. The 
theoretical results were corrected to agree with the wind-tunnel data in figure 1 in which 
y = 1.4 and Cj) = 1.61. The corrected Cj) curve is in reasonable agreement with the 
carbon tetrafluoride (CF4) data in reference 7 where for a y of about 1.16 to 1.17, Cj-j 
is in the range of 1.7 to 1.76. Since the flow behind the shock becomes supersonic when 
yg decreases to a value of approximately 1.12, theoretical Cj) values could not be 
obtained from the method of reference 6 below this value. Consequently, an arbitrary 
fairing from the Cj) value at y^ = 1.12 to the Newtonian value of Cj) = 1.766 at 
yg = 1 was employed. The use of perfect-gas results with yg from real-gas shock 
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FigTire 2.- Theoretical variation (plus empirical correction) of drag coefficient with 
effective ratio of specific heats for a 1^0® conical aeroshell. 

calculations is permissible since the ratio of specific heats is nearly constant in the shock 
layer ahead of extremely blunt bodies such as the 140^ cone. 

REAL-GAS THERMODYNAMIC PROPERTIES 


The free-stream real-gas effects are molecular rotation and vibration whereas for 
the higher temperature gas behind the shock molecular rotation, vibration, electronic 
excitation, dissociation, and ionization were considered. All equations and thermochemi- 
cal constants for computing thermodynamic properties from temperature, pressure, and 
elemental composition (relative amounts of C, O, and Ar of which the gas mixture is com- 
posed) are given in references 8 and 9. Surface conditions for the five model Mars atmo- 
spheres treated in the present paper ^Max Hp, Max Pg, Mean, Min Hp, and Min pg 
atmospheres^ are given in table I. 
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TABLE I.- ATMOSPHERIC CONDITIONS NEAR THE MEAN SURFACE FOR THE FIVE MODEL MARS ATMOSPHERES 





Conditions for — 



Parameter 







Min pg atmosphere 

Min Hp atmosphere 

Mean atmosphere 

Max Ps atmosphere 

Max Hp atmosphere 

Density scale height, 






km 

18,3 

11.9 

15.8 

10.8 

19.6 

ft 

60 000 

39 000 

51 800 

35 400 

64 300 

Density, 






kg/ra3 

0.00756 

0.0118 

0.0136 

0.0286 

0.0184 

g/cm3 

7.56 X 10-6 

1.18 X 10"6 

1.36 X 10-6 

2.86 X 10-6 

1.84 X 10-6 

slug/ft^ 

1.47 X 10-6 

2.29 X 10"6 

2.64 X 10-6 

5.55 X 10-6 

3.57 X 10-6 

Pressure, 






N/m^ 

400 

400 

600 

1000 

1000 

mb 

4 

4 

6 

10 

10 

Ib/ft2 

8.35 

8.35 

12.5 

20.9 

20.9 

Temperature, 

280 

180 

230 

180 

280 

Composition by mass, percent 






CO2 

100 

100 

88 

73 

73 

Ar 

0 

0 

12 

27 

27 

Composition by volume, percent 






CO2 

100 

100 

87 

71 

71 

Ar 

0 

0 

13 

29 

29 

Molecular mass 

44 

44 

43.5 

42.8 

42.8 

Acceleration of gravity, 






m/sec2 

3.72 

3.72 

3.72 

3.72 

3,72 

ft/sec2 

12,2 

12.2 

12.2 

12.2 

12.2 

Equatorial radius, 






km 

3393 

3393 

3393 

3393 

3393 

ft 

1.113 X 10'^ 

1,113 X 10'^ 

1.113 X 10'^ 

1.113 X 10'^ 

1.113 X lo'^ 


Free -Stream Thermodynamic Properties 

Free-stream temperature, pressure, and density are plotted against altitude in fig- 
ures 3 to 5 for the five models; static enthalpy and speed of sound were computed directly 
from free-stream composition (tables I and II) and temperature (fig. 3). The free-stream 
composition of each of the models is constant for altitudes up to 90 km (295 000 ft) and 
changes continuously at higher altitudes. Fortunately above 90 km, the density of four 
of the model atmospheres is so low that the variations in composition have negligible 
effects on the trajectories. The exception is the maximum density scale height model 
atmosphere (Max Hp atmosphere). Table II gives the altitude variation of composition 
for the Max Hp atmosphere, as table I gives only the low-altitude composition. 


Thermodynamic Properties Behind the Shock 

Polynomial approximations of the variations of equilibrium Ts and hg^RyT 3 
with pressure and density were used behind the shock wave. A given curve fit (poly- 
nomial approximation) is only valid for a constant elemental composition (constant propor- 
tions of the elements C, O, and Ar). As noted in the last section, the composition of each 
model atmosphere is constant for altitudes up to 90 km (295 000 ft), and only the Max Hp 
atmosphere is important at higher altitudes. 
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The Max Hp atmosphere required that equilibrium thermodynamic properties for 
three elemental compositions be used: one for the composition at altitudes up to 90 km 
(295 000 ft), another for the composition at 150 km (492 000 ft), and a third for that at 
240 km (787 000 ft). Since the elemental composition changes continuously above 90 km 
(295 000 ft), linear interpolation between the thermodynamic properties of these three 
altitudes was utilized to obtain the properties for the other altitudes. (See appendix C 
for details.) The interpolation parameter for hg/RyTig was the molecular mass of the 
undissociated gas mixture, whereas that for Tg was the number of elemental moles of 
carbon per unit mass of the atmosphere. Errors in hg^RuTg and Tg due to the inter- 
polation routine were usually less than 2 percent. 

RESULTS AND DISCUSSION 

Real-gas results are given for five model Mars atmospheres and a spectrum of 
entry conditions including orbital entry velocities of 4.27, 4.57, and 4.88 km/sec (14 000, 
15 000, and 16 000 ft/sec) and entry flight-path angles of 13. 5*^, 15°, and 17.5°. Entry 
conditions are defined at an arbitrary altitude of 244 km (800 000 ft). All calculations are 
for a 140° conical aeroshell with M/A = 76,5 kg/m2 (0.487 slug/ft2). Figures 6 to 14 
give calculated drag coefficients, dynamic pressures, and effective specific heat ratios 
as well as Mach numbers and flight -path angles. 

Drag Variation With Altitude 

The behavior of the real-gas drag coefficient as altitude changes is illustrated by 
figure 6(a) for the entry conditions Vg = 4,57 km/sec (15 000 ft/sec) and = 15°. 

For each model atmosphere, figure 6(b) locates the dynamic -pressure pulse or the alti- 
tude range where the value of the drag coefficient is most important. It is seen that at 
maximum dynamic pressure, the drag coefficients are in the neighborhood of Cp = 1.75 
which is about 9 percent higher than Cp = 1.61 in the wind-tunnel data of figure 1. 

The real-gas drag coefficient is, of course, determined from the variation of Cjj 
with shown in figure 2. Therefore, one could construct an altitude variation of yg 
from figures 6(a) and 2. For convenience, curves for the variation of yg with altitude 
are given in figure 7 and can be used together with figure 6(b) to show that yg is 
approximately 1.14 at maximum dynamic pressure. 

Drag Variation With Mach Number 

Figures 8 to 12 display the variation of drag coefficient with Mach number for all 
five model atmospheres and all entry conditions considered. The drag coefficient curves 
for the Max Pg, Mean, Min Pg, and Min Hp atmospheres all start with a value of 
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(b) I)ynamic pressure. 

Figure 6.« Variation with altitude of calculated drag coefficient and dynamic pressiore 
for = 15 *^ and V^, = ^*57 km/sec (15 000 ft/sec) for the five model atmospheres. 
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7.- Variation with altitude of effective ratio of specific heats for = 15° 
and Vjj = ^*57 km/sec (l5 000 ft/sec) for the five model atmospheres. 
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(a) Entry velocity effect for jZf-g = 15°. (b) Entry flight-path angle effect for 

Vg = 4.57 km/sec (15 000 ft/sec). 

Figure 8.- Variation of calculated drag coefficient with Mach number for the Max p model 
atmosphere for three entry velocities and three entry flight-path angles. 




(a) Entry velocity effect for = 15^. (b) Entry flight-path angle effect for 

V-g = 4.57 km/sec (15 000 ft/sec). 

Figure 9.- Variation of calculated drag coefficient with Mach number for the Mean model 
atmosphere for three entry velocities and three entry flight-path angles. 
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4. 88 16 000 

4.57 15 000 

4. 27 14 000 



(a) Entry velocity effect for ^ = 15 ^. (b) Entry flight-path angle effect for 

Vg = ^.57 km/sec (15 000 ft/sec). 

Figure 10 ,- Variation of calculated drag coefficient with Mach number for the Min p 
model atmosphere for three entry velocities and three entry flight-path angles . 



Cp. value for air, y - 1,4 



Cjj value for air. y = 1. 4 

\ I 


(a) Entry velocity effect for = 15 °, (b) Entry flight-path angle effect for 

Ve = ^-57 Wsec (15 000 ft/sec). 

Figure 11 .- Variation of calc-ulated drag coefficient with Mach number for the Min Hp 
model atmosphere for three entry velocities and three entry flight-path angles. 




(a) Entiy velocity effect for {2L = 15*^. (t>) Entry flight-path angle effect for 

Vg = 4.57 km/sec (15 000 ft/sec). 

Figure 12.- Variation of calculated drag coefficient with Mach numher for the Max Hp 
model atmosphere for three entiy velocities and three entry flight-path angles. 


approximately 1.78 at high Mach numbers and decrease monotonically to values varying 
from 1.62 to 1.67. Curves for the Max Hp atmosphere (fig. 12) are considerably dif- 
ferent in character from the others because drag coefficients for this model must be com- 
puted to much higher altitudes where the temperature and, hence, the speed of sound is 
high and the Mach numbers are low. Furthermore, the free -stream gas at high altitude 
is composed of atomic oxygen so that no dissociation takes place in the shock layer ahead 
of the body. Density jumps across the shock are therefore lower and this condition, in 
turn, leads to lower drag coefficients. The high-altitude ends of the Max Hp curves 
start at a Cj) of approximately 1.66 and a Mach number of about 8. As the altitude 
decreases, the drag coefficients and Mach numbers increase to a maximum and then 
decrease in a manner similar to the curves for the other atmospheres. 

It is evident from figures 8 to 12 that changes in entry velocity and flight-path angle 
over the ranges considered result in only small or negligible differences in the variations 
of drag coefficient with Mach number. The largest effect of entry velocity changes occurs 
in the drag coefficient curves for the Max Hp atmosphere whereas small increments in 
the flight-path angle at entry have the most effect in the Min Hp atmosphere. 

Drag coefficients for all five model atmospheres are compared in figure 13 for the 
entry conditions Vg = 4.57 km/sec (15 000 ft/sec) and = 15^. Ideally, one would 


17 




Figure 15-- Variation of calculated drag coefficient with 
Mach number for an entry flight-path angle of 15° and 
an entry velocity of 4.57 km/sec (15 000 ft/sec) for 
all five model atmospheres. 

hope to find a single Mach number variation of drag coefficient that would be satisfactory 
for all Mars orbital entry trajectory calculations. If such a curve did exist, it would 
afford a substantial convenience in parametric studies. The magnitude of error that one 
would incur by using an average drag curve can be obtained from figure 13 which shows 
that drag coefficients increase approximately 3 percent in going from the Max pg to the 
Min Pg atmospheres. An average real-gas drag curve will result then in roughly a 
1.5-percent error. In the light of the uncertainties in the Mars atmosphere itself, the 
selection of a single drag coefficient curve, be it an average or a conservative one, for 
mission analysis work may be desirable. 

Trajectory Comparisons 

To illustrate the difference between real-gas drag coefficient trajectories and those 
obtained by using a constant drag coefficient of Cj) = 1.61, figure 14 and table HI have 
been prepared. Figure 14 shows the variations of Mach number and flight-path angle with 
altitude for the five model atmospheres. Inasmuch as some assumptions inherent in the 
present analysis become less valid with decreasing Mach number, the computations were 
terminated at a Mach number of 2. 

Since one would expect any significant differences in the trajectory parameters to 
appear most pronounced after integrating over the entire trajectory (down to Mach 2 in 
the present calculations), pertinent numerical results are given in table in. It is seen 
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15 f^eg 
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Altitude 


(a) Max pg model atmosphere. 



(b) Mean model atmosphere. 

Figiore l4,- Altitude variation of Mach number and flight -path angle for entry- 
conditions of Vg = ^.57 km/sec (15 000 ft/sec) and - 15 ° for the 

five model atmospheres. (See table HI for a tabulation of trajectory- 
parameters at a Mach number of 2 . ) 
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M 
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8h 
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(e) Max Hp model atmosphere. 
Figure l4.- Concluded. 


TABLE m.- REAL-GAS AND CONSTANT -DRAG COEFFICIENT DETERMINATIONS OF V^, 0, AND H AT MACH 2.0 

[Ve = 4.57 km/sec (15 000 ft/sec); 


Model 

atmosphere 

Constant drag 

v„ 

Real -gas drag 
Vco 

Constant 
drag 0, 
deg 

Real -gas 
drag 0, 
deg 

Constant drag 
H 

Real -gas drag 
H 

Real-gas drag H 
minus constant 
drag H 

m/sec 

ft /sec 

m/sec 

ft/sec 

m 

ft 

m 

ft 

m 

ft 

Max Pg 

399.0 

1309 

398.4 

1307 

21.79 

22.15 

14 910 

48 910 

15 090 

49 500 

180 

590 

Mean 

456.9 

1499 

456.3 

1497 

23.00 

23.30 

13 430 

44 070 

13 850 

45 450 

420 

1380 

Min pg 

489.8 

1607 

490.4 

1609 

23.81 

24.01 

10 700 

35 090 

11 340 

37 210 

650 

2120 

Min Hp 

391.1 

1283 

389.5 

1278 

21.94 

22.25 

7 379 

24 210 

7 654 

25 110 

270 

900 

Max Hp 

506.9 

1663 

506.3 

1661 

23.26 

23.58 

22 900 

75 130 

23 370 

76 680 

470 

1550 
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therein that there are significant differences in altitudes but very small differences in 
flight -path angles and velocities. The altitude differences range from 180 meters (590 ft) 
for the Max Pg atmosphere to 650 meters (2120 ft) for the Min pg atmosphere. For 
the Min Hp atmosphere, where a Mach number of 2 is reached at the lowest altitude, 
real-gas effects provided a 270-meter (900-ft) increase in altitude. 

In connection with the Max Hp trajectory in figure 14(e) it should be noted that a 
substantial effort (see appendix C) was expended on the real-gas calculations above 90 km 
(295 000 ft). Subsequent simplified calculations using a constant Cj^ of 1.78 above 
90 km (295 000 ft) and the real-gas Cj) at lower altitudes yielded a trajectory which 
was in excellent agreement with the real-gas Max Hp model trajectory. It appears 
that the application of this simplified approach in determining trajectories for model 
atmospheres similar to the Max Hp atmosphere would be justified. 

CONCLUDING REMARKS 


The real-gas drag coefficient was found to vary from about 1.62 to 1.78 and to have 
a value near 1.75 at a maximum dynamic pressure for orbital entry into the five model 
Mars atmospheres considered; 


Max Hp atmosphere 
Max Pg atmosphere 
Mean atmosphere 
Min Hp atmosphere 
Min Pg atmosphere 


atmosphere with maximum 
atmosphere with maximum 
mean model atmosphere 
atmosphere with minimum 
atmosphere with minimum 


average density scale height 
density near surface 

average density scale height 
density near surface 


For a given model atmosphere, very small differences in the drag coefficient occurred 
when entry conditions were varied over the ranges considered. However, the difference 
between the Mach number variation of drag coefficients for the Min pg and Max pg 
model atmospheres is about 3 percent and the other three generally fall between these 
two. The real-gas variation of drag coefficient with Mach number for the Max Hp 
model atmosphere was different in character from the variations of the other four model 
atmospheres. 

The differences between real-gas drag coefficient trajectories and those obtained 
by using a constant drag coefficient of 1.61 were shown. At a Mach number of 2, there 
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were very small differences in flight-path angles and velocities but significant differences 
in altitudes, the Min pg model atmosphere giving the largest difference (650 meters or 
2120 ft). A Mach number of 2 is reached at the lowest altitude for the Min Hp model 
atmosphere, where real-gas effects provide a 270-meter (900-ft) increase in altitude. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., January 21, 1971. 
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APPENDIX A 


EQUATIONS OF MOTION 


For completeness, the standard point mass trajectory equations for a nonlifting 
body are given in this appendix. The derivatives of the velocity, flight-path angle, and 
altitude are 


dV. 

dt 


= g sin <p 




M/CpA 


(Al) 


d(p 

dt 


1 ( yJ 

^ocVS + H 


glcos 


(A2) 


dt 


-V^ sin (p 


(A3) 


where 


APPENDIX B 


SHOCK DENSITY RATIO DETERMINATION 


The purpose of this appendix is to show how the shock density ratio Pq/Poo 
computed from the instantaneous velocity and ambient conditions in the free -stream and 
equilibrium enthalpy as a function of pressure and density behind the shock. The 
oblique-shock equations (ref. 5) are solved by a procedure similar to that of reference 10. 

The component of the free-stream velocity normal to the shock wave is 

Uoo = sin e (Bl) 

where 6 is a prescribed shock inclination angle. An initial estimate is made for the 
density behind the shock pg and is used to compute the normal component of the veloc- 
ity behind the shock: 

Poo" 00 (B2) 

S Pg 

The pressure and enthalpy behind the shock based on the estimated value for pg are 
computed by using equations (Bl) and (B2) and the ambient conditions: 

PS = + Poo + Poo"oo^ 




oo^ ^oo 


(B4) 


Equilibrium enthalpy behind the shock is evaluated from polynomial curve fits of its 
variations with pressure Pg and density Pg. Since Pg is an estimated quantity, the 
enthalpy thus obtained agrees only approximately with that given by equation (B4). Addi- 
tional values of density are substituted into the thermodynamic curve fit equations until 
the value of hg obtained agrees with equation (B4) to within 100 mVsec? (1076 ft^/sec^). 
The value of density which provides the aforementioned enthalpy agreement is then sub- 
stituted back into equation (B2) and the process is repeated. When two successive values 
of Pg (that is, those required to obtain enthalpy agreement) agree to within 10"^ per- 
cent, the calculation is terminated. 
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APPENDIX C 


fflGH-ALTITUDE THERMODYNAMIC PROPERTIES FOR 
MAXIMUM Hp ATMOSPHERE 

The interpolation routine which computes the equilibrium properties hg^RyTs 
and Tg at altitudes from 90 to 240 kilometers (295 000 to 787 000 ft) for the Max Hp 
model atmosphere is described in this appendix. Three sets of curve fits (polynomial 
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4 . 200 
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Figure 15*- Variation of temperature with the interpolation parameter ^^C02 
densities at Pg/Po ~ 10“5 for the Max Hp atmosphere. 


26 



APPENDIX C - Concluded 


approximations of thermodynamic properties) are used, valid for the atmospheric com- 
positions at 90, 150, and 240 kilometers (295 000, 492 000, and 787 000 ft). For a given 
pressure ratio and density ratio, each property is determined by linear interpolation from 
two of the three sets of thermodynamic -property curve fits. The interpolation parameter 
for hs^RyTg is the molecular weight My of the undissociated gas mixture, that is, 
the mixture which would exist at low temperatures where all the CO and O are combined 
into CO 2 and O 2 . The interpolation parameter for Tg is Nco2 which is the number 
of moles of the element carbon per unit mass of the atmosphere (equal to the number of 
moles of CO 2 per unit mass of the undissociated mixture). These parameters were com- 
puted from the composition of the Max Hp model atmosphere (table n) and tabulated at 
10-kilometer (32 800-ft) intervals for altitudes of 90 kilometers (295 000 ft) and higher. 

Figures 15 and 16 show how well the linear interpolation procedure works for a pres- 
sure ratio of Pg/Po = lO"^ and two values of density ratio, Pg/P^ = 10“^'^ and 10“^*^. 
The circled points are the exact values (ref. 8) whereas the straight-line segments repre- 
sent the interpolation values. Circles were not placed around the points at altitudes of 90, 
150, and 240 kilometers (295 000, 492 000, and 787 000 ft), since the straight lines were 
drawn by connecting these points. The differences between the exact values (circles) 
and the interpolated values (straight lines) are usually less than 2 percent. 



Figure 16.- Variation of the thermodynamic quantity hg/RuTg with the 
interpolation parameter for two densities at Ps/P© “ 10“^ for 

the Max Hp atmosphere. 
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